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Due to its good economy and applicability, high-power thyristor device is widely used in power supply sys- 
tems of large fusion devices. The high-dose neutrons produced by Deuterium-Tritium (D-T) fusion reactions 
will irradiate the thyristor for a long time, which will lead to the change of the electrical characteristics of the 
thyristor, and may eventually cause irreversible damage. In this paper, the thyristor switch of the commutation 
circuit in Quench Protection System (QPS) of fusion device is regard as the study object. The relationship be- 
tween the internal physical structure and the external electrical parameters of irradiated thyristor is established. 
Then, a series of targeted thyristor physical simulations and neutron irradiation experiments are carried out to 
verify the correctness of the theoretical analysis. In addition, the effect of irradiated thyristor electrical charac- 
teristic changes on the whole QPS is studied by accurate simulation, which provides valuable guidance for the 


maintenance and renovation of QPS. 


Keywords: Fusion device, Neutron Irradiation effects, Thyrisor, Quench protection 


I. INTRODUCTION 


Tokamak is an experimental device for the development of 
nuclear fusion energy, in which the most easily achieved fu- 
sion reaction is Deuterium-Tritium (D-T) fusion, whose reac- 
tion formula is shown in Eq. (1) [1, 2]. When D-T fusion 
reaction occurs in Tokamak, about 80% of the energy will be 
released in the form of 14.1 MeV fast neutrons [3-5]. Nowa- 
days, Tokamak is still in the experimental exploration stage, 
and some windows and channels need to be set up for obser- 
vation, which may lead to an increase in neutron irradiation 
flux in fusion environment. Taking the ITER Tokamak as an 
example [6], when the 500 MW D-T fusion reaction experi- 
ment is carried out, the number of fast neutrons produced by 
the plasma is 1.78 x 107°n x s-t [7]. After installing the 
shielding layer and considering the neutron attenuation, the 
calculated neutron flux rate in the whole host hall is still up 
to 10°n x cm7? x s71, which is far great than the safety 
threshold and has drastic impacts on the performance of the 
electrical devices around Tokamak [8, 9]. 


D +T > n(14.1MeV) +“He (3.5MeV) , (1) 


Quench Protection System (QPS) is usually designed to en- 
sure the safety of the superconducting magnet in Tokamak 
[10, 11], in which thyristors are commonly used as the trigger 
unit of the commutation circuit, as shown in Fig. 1 [12, 13]. 
QPS works as follows: when superconducting magnet in nor- 
mal operation condition, the current flows through ByPass 
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27 Switch (BPS). Once quench occurs, BPS quickly opens, and 
23 the magnet current is transferred to Vacuum Circuit Breaker 
29 (VCB) by the arc voltage. Then the thyristors turn on, a pulse 
30 Current generated by charged capacitor C and pulsed inductor 
31 L, to force the current in VCB down to zero [14, 15]. Finally, 
32 the magnet current is transferred to Discharge Resistor (DR) 
33 and the energy stored in magnet is constantly consumed until 
34 to zero. 


Discharge Resistor 


Commutation 
Circuit 


Fig. 1. (Color online) The circuit topology of QPS. 
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in the fusion device, the QPS needs to respond as quickly as 
possible. To ensure the QPS can do this, it should be located 
as closely as possible to the magnetic coils. Unfortunately, 
this is also where there is a high fast neutron flux rate, thus 
we likely cannot avoid neutron damage. The typical layout of 
fusion device and irradiation nephogram are shown in Fig. 2 
[8, 9]. If the thyristor switch is damaged due to neutron irra- 
diation and cannot operate normally, the current would not be 
interrupt quickly by VCB with the help of the auxiliary pulse 
current, which resulting in the failure of QPS and the destroy 
of superconducting magnet [16]. 
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Fig. 2. (Color online) Typical layout of fusion device and irradiation 

nephogram. (a) Typical layout of fusion device. (b) Total Neutron 

Flux (TNF) caused by plasma in host hall, the max TNF is 10°n x 
-2 “1 

cm xS. 


Since ITER takes the D-T fusion route, neutron irradiation 
damage is an issue that must be considered. More previous 
studies are focused on the main structural materials of fusion 
reactor, such as the First Wall, Blanket Modules and Divertor 
[17-19]. Gradually, ITER realized the threat of neutron irra- 
diation to auxiliary equipment, and conducted irradiation cal- 
culation, established exposure limits and provided mitigation 
measures [8, 20, 21]. Researchers have carried out radiation 
reliability analysis on the basis of ITER guidelines, but it is 
far from sufficient, especially for high-power thyristors. 
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Semiconductor damage via high energy neutron irradiation 
is a topic not only of interest in the context of Tokamaks and 
Stellarators, but also in the context of spacecraft and satel- 
lites, where the electronics are similarly bombarded with a 
complex mix of high energy particles. [22—24]. Currently, in 
literature, the displacement damage effect and total ionization 
effect occurred in electronics and few power device have been 
separately studied at macro and micro levels. However, as the 
particularity of fusion neutron irradiation environment and 
the requirement for effective irradiation reinforcement, the re- 
lationship between internal physical structure changes and ex- 
ternal electrical characteristic changes should be deeply stud- 
ied and accurately built. 

The neutron irradiation damage and internal physical struc- 
ture change of semiconductor devices are mainly carried out 
by experiments. S. Yue et al. studied the synergistic effect 
of electrical stress and neutron irradiation on silicon carbide 
power MOSFETs. The defects in the devices were charac- 
terized by the Low Frequency Noise (LFN) and Deep-Level 
Transient Spectrum (DLTS) methods [25]. R. Chen et al. car- 
ried out the Displacement Damage Dose (DDD) and Total 
Ionization Dose (TID) irradiation experiments with the 14.1 
MeV neutron high-voltage multiplier and 60Co gamma-ray 
to study the synergistic effect of TID and DDD in neutron ir- 
radiated GaN HEMT [26]. However, the change relationships 
between various physical parameters and electrical character- 
istics of semiconductor devices need to be further studied. 


On the other hand, F. Ravotti and J. Mekki studied I- 
V characteristic of PIN photodiode under 1 MeV neutron 
irradiation at European Organization for Nuclear Research 
(CERN) [27, 28]. Through a lot of irradiation experiments 
and tests, the relationship curves of forward conduction volt- 
age drop and leakage current under different neutron flux are 
obtained. C. Liu ef al. studied the influence of fast neu- 
tron irradiation on the high di/dt switching characteristics of 
the pulse power supply Insulated Gate Triggered Thyristor 
(IGTT) through experiments [29]. The simulation and test 
results indicate that the neutron irradiation bring significant 
decrease of threshold voltage and on-state voltage of IGTT. 
However, the studies did not give a detailed description and 
exploration of the physical causes leading to the electrical 
characteristic changes. 


In this paper, we focus on studying the characteristics 
change of the thyristor switch within the fusion device’s QPS 
shown in Fig. 1. The internal physical structure changes of the 
thyristor caused by neutron irradiation are introduced firstly 
to obtain the possible physical parameters which may cause 
the electrical characteristic changes. Then, in view of vari- 
ous key electrical parameters of thyristor, the relationship be- 
tween physical and electrical characteristics of thyristor after 
neutron irradiation is theoretically analyzed. Also, a series 
of targeted thyristor physical simulations and neutron irradi- 
ation experiments are carried out to verify the correctness of 
the theoretical analysis results. Furthermore, the effect of ir- 
radiated thyristor characteristic changes on the whole QPS is 
also been studied by the performance simulation. The study 
not only build a clear relationship between the changes of 
physical and electrical characteristics of the neutron irradi- 
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ated thyristor, but also provides a valuable guideline for the 
maintenance and renovation of QPS. 


II. THEORETICAL DISCUSSION OF THE 
RELATIONSHIP BETWEEN PHYSICAL AND 
ELECTRICAL CHARACTERISTICS CHANGE OF THE 
THYRISTOR AFTER NEUTRON IRRADIATION 


A. Damage Mechanism Analysis of Physical Characteristics 
of Thyristor after Neutron Irradiation 


Since neutrons are electrically neutral and have very strong 
penetrability, it can be sufficiently close to the atomic nucleus 
of the lattice atom of the material to be radiated and collide 
with the atomic nucleus elastically [30, 31]. The lattice atom 
gets energy in the collision process, thus leaving its normal 
lattice position and becoming an interstitial atom in the lat- 
tice, and leaving a vacancy in its original position. If the va- 
cancy and interstitial atom are still within the coulomb force 
field of its elastic force field, recombination can be occurred. 
If the force field is exceeded, the interstitial atoms will not 
recover, thus causing atomic displacement. The production 
process of vacancy and interstitial atom is shown in Fig. 3. 


Incident 
Neutrons 


Fig. 3. (Color online) The production mechanism of vacancy and 
interstitial atom in Si material. 


Vacancies can be combined with adjacent atoms or with 
vacancies. Vacancies can also be moved to the vicinity of 
doping to form vacancy impurity complexes, so that doping 
do not participate in conductivity, thus changing the density 
of doping. The displacement of a single atom forms a simple 
defect, called the Frenkel Defect [32]. If the incident neutrons 
energy is large enough, the displaced interstitial atoms will 
get enough energy. The collision of such interstitial atoms 
will displace a large number of atoms in the lattice, forming 
a large defect cluster. 

Literature [33] indicates that when the incident neutron en- 
ergy is 1 MeV, the average energy obtained by the initially 
displaced atom is 72.5 keV, while the displacement threshold 
for a Si atom is 15 eV. The D-T reaction produces a neutron 
energy of 14.1 MeV, which is about 1 million times the Si 
atomic shift threshold, thus obviously has a great impact on 
Si materials. In addition, it also shows that a defect group 
generated by a 1 MeV neutron can include about 200 to 1000 
displaced lattice atoms, Although the collision cross section 
between neutron and lattice atom is only between 10723 cm? 
and 10~*cm?, the energy transferred to the lattice atom by 
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collision is quite large, a series of lattice atoms can eventually 
be displaced. 


Semiconductor materials have three microscopic physical 
parameters, namely minority carrier lifetime 7, pure doping 
density n and mobility u. Neutron irradiation introduces de- 
fects and defect cluster in the semiconductor, thus changing 
the three microscopic physical parameters and changing al- 
most all parameters of the semiconductor device. 


Defects formed by neutron irradiation can introduce addi- 


65 tional energy levels into the forbidden band of semiconduc- 


tor materials, which increases the chance of recombination of 
minority carriers and majority carriers, increases the recom- 
bination rate, and thus reduces the minority carrier lifetime. 
The relationship between neutron flux and the minority car- 
rier lifetime is given below. 


oes 


T (9) 5 


Where 7 (0) and 7 (@) are the minority carrier lifetime be- 
fore and after neutron irradiation, k, is the minority carrier 
lifetime damage coefficient of semiconductor materials, is 
the neutron flux. 


As thyristor is a kind of semiconductor device with minor- 
ity carriers as conductive mechanism, thus minority carrier 
lifetime damage will greatly affect the performance of the de- 
vice. 

Beside this, neutron irradiation also reduces the doping 
density of the semiconductor materials, which called carrier 
removal effect. The relationship between neutron flux and the 
doping density is given below. 


(3) 


Where n (0) and n(@) are the doping density before and 
after neutron irradiation, An/A¢@ is the carrier removal rate, 
and ¢ is the neutron flux. 

Also, when neutron irradiation introduces defects in semi- 
conductor materials, these defects can be used as the scatter- 
ing center of carriers, effectively reducing the mobility of car- 
riers. The influence of these neutron irradiation on mobility 
can be expressed as 


1 A, x [1 —exp (er — es) /kT] 
+ 
HL Bi 


d, 
(4) 


Where u (@) is the mobility of carriers after neutron irra- 
diation, urp is the mobility due to stray atom scattering, wr 
is the mobility due to collision with thermal vibration lattice, 
A, and B; are the coefficients, €r and € are the defect level 
and Fermi level respectively, k is Boltzmann constant, T is 
thermodynamic temperature, and ¢ is the neutron flux. 


200 


B. Establishment of Relationship between Physical and 


21 Electrical Characteristics Change after Neutron Irradiation 


22 As thyristor has longer lifetime, more reliability, less main- 
203 tenance and high rated power capability compared to other 
204 Semiconductor devices, so it becomes the first choice for QPS 
205 in Tokamak superconducting system. 
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Fig. 4. (Color online) Internal structure diagram of thyristor. (a) 
Blocking state. (b) Conducting state. 


Thyristor is composed of P-type and N-type semiconductor 
alternately [34]. As shown in Fig. 4, the depletion layer is 
formed at the junction of P-type and N-type semiconductor. If 
there is no external power supply to the gate of thyristor, there 
will always be at least one reverse biased depletion layer re- 
gardless of the voltage direction of the external circuit, which 
are shown in Fig. 4(a). If the thyristor gate is applied by an 
additional voltage, which is called gate trigger, a large num- 
ber of electrons will be injected into P-region of the thyristor, 
and P-region eventually become N-region which will cause 
the depletion layer width gradually decrease and eventually 
almost disappear as shown in Fig. 4(b). As the injected elec- 
trons have already entered N-region, even though the thyristor 
gate is not applying additional voltage, it can still maintain the 
conducting state. 
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Hence, both blocking and conducting state of thyristor after 
neutron irradiation can be analyzed by simplifying it to PN 
diode. 


221 


222 


223 


224 1. Blocking Characteristic 


25 The blocking charactersitic of thyristor is actually an ex- 
226 pression of the reverse current of PN junction under reverse 
227 voltage, which can be divided into two conditions given in 
28 Fig. 4(a). Hence, the forward leakage current prm and the 
229 reverse leakage current I Rrm are the most two crucial param- 
230 eters to evaluate the blocking characteristic of thyristor. Fig. 
231 5 shows the change of PN junction barrier region with reverse 
232 bias. When the applied electric field E is in the same direc- 
233 tion as the self-built electric field Eo, the electric field in the 
234 depletion layer is strengthened. The enhancement of electric 
235 field makes the barrier region change in two ways: first, the 
23 barrier region becomes wider, and second, the barrier height 
237 increases. 


Ent 


Fig. 5. (Color online) The change of PN junction barrier region with 
reverse bias. 


La 


Fig. 6. (Color online) The diagram of reverse current generation and 
minority distribution. 


28 Meanwhile, the balance between diffusion and drift of the 
233 carriers is broken by the growth of the barrier region, and drift 
240 plays a dominant role. Therefore, the electrons at the bound- 
24 ary of the P-region could be pulled to the N-region, and the 
242 holes at the boundary of the N-region could be pulled to the 
243 P-region. In the meantime, the electrons in the P-region and 
244 the holes in the N-region will run to the boundary to supple- 
24s ment, thus forming the reverse current J, which is directed 


24 from the N-region to the P-region, as shown in Fig. 6. How- 
247 ever, Since there are few electrons in the P-region and holes 
248 in the N-region, the reverse current Ip is very small. 

29 In the case of ignoring the recombination of electrons and 
250 holes in the barrier region, the current density J through the 
21 PN junction can be considered as the sum of the hole dif- 
2s2 fusion current density J, (~~) passing through the zy, the 
23 electron diffusion current density J, (xp) passing through 
the x p [34]. 


254 


255 J = Jp (tn) + Jn (zp), (5) 


26 J, (xn) could be calculated by the following equations if 


257 the hole and electron density distribution are obtained. 


dApn (x) 
dx , 


L=IN 


Jy (tn) = —qDp (6) 
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29 Where dAp/dz are the hole density gradient, D, is the 
hole diffusion coefficient. 

If the P and N regions are long enough, the non-equilibrium 
carrier concentration in the PN junction should be distributed 


exponentially. 
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Apn (£) = [Pn (an) — Pro] exp (= — z) E 


Where Lp is the hole diffusion length. 
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266 Lp =v DpTp, (8) 
267 Also 

U 
200 Pn (EN) = Pno XP (5) (9) 


Where pno is the equilibrium minority carrier concentra- 
tion in N-region, U is the applied reverse voltage. Hence, 
Jp (xy) could be calculated by the following equation. 


qDpPno qU 
272 Jp (zN) = a lex» (5) = 1 5 (10) 
273 By the same token, 

Dn U 
253 Consequently, 

qDnnpo , qDpPno qU 
= —— ]-1 12 

276 J ( L, + o ) lex (5 ’ (12) 


qDnnpo 


a 


qDpPno qU 
Ip = AT = A( ) Jew (2) il 
(13) 
23 Where A is PN junction cross-sectional area. As U ap- 
279 proaches negative infinity, exp (qU /kT) approaches 0. 


qDnNpo 
A 
Z 
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qDpPno 


I 
F in 


+ 


280 ) = —Ip, (14) 

Hence, the J would increase as the carriers’ lifetime Tp 
and 7,, would be greatly shortened after neutron irradiation. 
In addition, as for silicon PN junction, the recombination cur- 
rent density J; in barrier region also plays a crucial role in 
the reverse current. After neutron irradiation, the number of 
recombination centers in the barrier region of PN junction in- 
creases, and the carrier generation rate increases, which leads 
to the rise of recombination current. 

As a result, after neutron irradiation the reverse current Ip 
in PN junction would be increased, which will result in the in- 
crease of forward leakage current Ip ray and reverse leakage 
current [pry of irradiated thyristor. 
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293 2. Conducting characteristic 


The conducing characteristic of semiconductor can be 
characterized by the on-state resistance or on-state voltage 
drop Urm. The lower the Urm, the better the high current 
performance of the device, which is a significant parameter of 
the semiconductor applied in QPS of fusion device as the ex- 
tremely high current breaking requirement of superconduct- 
ing magnet. 

In order to have high working voltage capability, the sub- 
strate resistivity of the PN junction diode is generally high, 
that is, the substrate doping concentration is low and the sub- 
strate region is wide. Therefore, the defects caused by neu- 
tron irradiation led to carrier removal effect in semiconductor 
material is more apparent, resulting in plenty of conductive 
catrier constrain and not participate in conduction. Hence, 
the resistivity increases, so as to increase the conducting volt- 
age of PN junction. 

According to the above analysis, the conducting charac- 
teristic of thyristor can be analyzed by simplifying it to PN 
diode. The on-state voltage drop of PN diode Urm can be 
divided into two parts: one is the junction voltage U; formed 
by the internal electric field, and the other is the voltage re- 
lated to the substrate epitaxial layer Up, U; and is very small 
and negligible. 


Ur =U; + Up ~ Ub, (15) 


317 


sis According to the literature [33], Up can be expressed as the 
319 following equation. 


Iwa 


[p= u 
b Aon’ 


(16) 


set Where wg is the base region width of PN junction, op is 
the conductivity of P-region. 


The conductivity of irradiated material can be derived as 


322 


323 


a (¢) =a (0)e"°?, (17) 


324 


o (0) 
325 In = ad, (18) 
a (9) 
3 When ¢ approaches 0, 
1 
2 v= Te), (19) 
328 As 
329 Oo = qHP, (20) 


Where u is majority carrier mobility of PN diode, and p is 
base doping concentration of PN diode. 


332 Then 
333 R= X ; (21) 
do o—0 
334 a= ue x R, (22) 
o(0) _ qu(d) 
335 = = x Rọ, (23) 
o(¢) (0) 
336 — Where R is carrier removal rate. 
37 According to Eq. (16) 
I 
a o= o (24) 
39 Substitute Eq. (24) into Eq. (23) 
Up () qu ($) 
= 1 2 
340 U, (0) exp o (0) x Rọ| > 1, (25) 


31 It can be concluded from the above analysis that U, will 
342 increase as majority carrier removal effect in the PN junc- 
343 tion. In addition, the change of the irradiated thyristor on- 
344 State voltage drop Urm is determined by the irradiation in- 
345 tensity, which will be discussed by further neutron irradiation 
34 Simulation and experiment. 


347 3. Reverse Recovery Characteristic 


34a Reverse recovery is “turning off’ a thyristor after it has 
34 been in a conducting state, and that the time taken for it to 
350 SWitch from conducting to blocking state is critical and needs 
351 to be minimised. Usually neutron damage decreases reverse 
352 recovery time, this is because the minority carrier lifetime of 
353 the thyristor is significantly shortened after irradiation. 

34 The typical reverse recovery current waveform is shown in 
355 Fig. 7. Ip is the peak value of the forward current. ts = 
356 tg — tı is the time of the reverse current grows from zero to 
37 Maximum value Irm, which is called storage time. ty 
358 t3 — tə is the time of the current drop from Ipyy to O.L Ry, 
359 Which is called fall time. The total reverse recovery time is 
360 trr = ts — ty, and Qrr is the reverse recovery charge [35]. 


Fig. 7. (Color online) The typical reverse recovery current wave- 
form. 


31 When the thyristor is in conducting state, most of the 
32 Charges stored in the base region, which are formed by the 
363 Minority carrier. The charge continuity equation is given as 
364 follow, 


_ dQ, Q 


ir (t) = TE (26) 


365 


Where ir (t) is the current flow through the thyristor, Q 
is the stored charge in base region, 7 is the minority carrier 
lifetime. 

The thyristors in QPS play the role of the trigger unit of the 
commutation circuit composed of the pulsed inductor L and 
the charged capacitor C. To simplify the analysis, the for- 
ward conducting current is considered to be sinusoidal cur- 
rent, ip (t) = Ipsinwt. According to Eq. (26), it can be 
obtained that 


366 
367 
368 
369 
370 
371 
372 
373 


374 


Kol 
Q(t) = eee (wre + sin wt — wT cos wt) ; 


(27) 


375 
376 


Where Ko is the current gain of the thyristor. 
At tı, the charge Q (t1) can be calculated by using Eq. (27) 


377 


o Kolpwr? 


Q (ti) = EE (28) 


378 


(1 + ens) : 


379 | When the forward current drops to zero, there are still large 
number of non-equilibrium carriers in thyristor, which will be 
dissipated by the way of sweeping out and recombination. At 
time t2, the current reaches Ipm, and the thyristor begins to 
recover the reverse blocking capability. At this time, only a 
small amount of stored charge remains in the base region. 
Those charge will disappear by the way of recombination, 
which is classified as low injection condition. Actually, high- 
387 power semiconductor devices generally work under large in- 
338 jection condition, so it can be approximately considered that 
339 Q (t2) = 0. According to Eq. (27), it can be calculated that 


380 
381 
382 
383 
384 
385 
386 


390 sin wtp = wT (cos wt — gta) ; (29) 


31 As t2 = tı +t, = T/w + ts, wte = T + wts. Also, wts is 
32 far less than 7/2, so wts approximately considers to be zero. 
303 Then, it could be deduced that sin wt, ~ wts, coswt, œ 1, 
394 therefore, 


395 sin wtp = sin (t + wt,) = —sin(wt,) ~ —wt,, (30) 


396 cos wta = cos (m + wts) = — cos (wts) ~ —1, (31) 
Approximating Eq. (29), the expression of tə can be ob- 


tained 


397 


398 


399 (32) 


ts ST (1 +e¥/7) ; 


40 t2 is much longer than 7, so 


401 ts © T, (33) 
Then, the maximum reverse current Irm can be expressed 


by the following equation 


Irm ~ k (to —t1) = kz, (34) 


404 


45. Where k is the forward current change rate at tı point. 

On the basis of above analysis, it can be concluded that af- 
ter high-dose neutron irradiating the thyristor, which reduces 
the minority carrier lifetime 7, the reverse recovery charac- 
teristics will have significant changes. According to Eq. (33) 
and Eq. (34), the stored time t, and the maximum reverse cur- 
rent Irm will decreases with the reduction of 7, which means 
that the total turn-off time of the semiconductor tyy will be 
shorten as the result of high-dose neutron irradiation. Also, 
the reverse recovery charge Q,,. will have a sharp decrease 
after irradiation. 
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II. NEUTRON IRRADIATION SIMULATION 


416 


“7 To verify the correctness of the above theoretical analysis 
of the thyristor damage mechanism under fusion neutron irra- 
diation, a qualitative verification through material and device 
level simulation based on simplified thyristor model is neces- 
sary. The idea of simulation in this paper is that after neutron 
irradiation, the lattice of the device material is damaged, and 
the defects introduce additional energy levels, which will af- 
fect the electrical characteristics of the device. Therefore, two 
types of defects are mainly considered in this paper, namely, 
donor defects, characterized by Æc-0.12 eV, and acceptor de- 
fect, characteristic Ey +0.21 eV. 

APSYS, as an advanced semiconductor device physical 
model simulation software, is used for simulation in this pa- 
per. The thyristor is an extremely complex device, and its 
designs are often patented and proprietary. Hence, these de- 
signs are guarded commercial secrets. Therefore, a general 
and widely used physical model of thyristor is established 
to simulate the effect of neutron irradiation on the electrical 
characteristics of thyristor devices, the 2D model of the de- 
vice and doping concentration as shown in Fig. 8. 
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Fig. 8. (Color online) Simulation model of thyristor. (a) 2D model. 
(b) doping concentration. 


The mathematical models mainly include basic models and 
ass key models. The basic model include drift diffusion model, 
poisson equation, continuity equation, ohmic contact electri- 
cal boundary, etc.The key models include mobility model, 
collision ionization model, bulk defect model, etc. After de- 
termining the physical parameters of the thyristor model, the 
basic semiconductor equations can be solved numerically by 
combing the initial conditions and the boundary conditions, 
and the electrical parameters under the specific bias voltage 
are solved. 

The blocking characteristic of the thyristor are shown in 
Fig. 9a. It can be seen that the leakage current of the thyristor 
increases sharply when the reverse voltage is 350 V. When 
Uax = -400 V, the 2D electric field distribution is shown in 
Fig. 9b, the main electric field falls in the N-region. 

The on-state characteristic of the thyristor is shown in 
453 Fig. 10. It can be seen that when the triggering voltage Vg 
454 = 0.5 V, the thyristor does not trigger on, while when Vg = 1 
455 to 5 V, the thyristor triggers on. 

The comparison of the blocking characteristics of thyristor 
457 before and after neutron irradiation are shown in Fig. 11 and 
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Fig. 9. (Color online) Blocking characteristic. (a) I-V characteristic. 
(b) 2-D electric field distribution (U a x = -400 V). 


Fig. 13. (Color online) Urm of irradiated thyristor (donor-type de- 
fect). (a) Linear. (b) Semi-log. 
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Fig. 10. (Color online) Conducting characteristic. (a) Linear. (b) 
Semi-log. 
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Fig. 14. (Color online) Urm of irradiated thyristor (acceptor-type 
defect). (a) Linear. (b) Semi-log. 


The comparison of on-state characteristics of thyristors be- 
fore and after neutron irradiation is shown in Fig. 13 and 
Fig. 14. It can be seen that when the thyristor flows through 
the same current, the on-state voltage drop Urm increases 
after irradiation under both donor defect and acceptor defect 


Semi-log J 
100 200 300 aoo 47 effects. l , 
Vax (V) Vax (V) ass The simulation results show that the forward and reverse 
(a) (b) 49 leakage current and on-state voltage of thyristor increase after 


Fig. 11. (Color online) Jp rm of irradiated thyristor (donor-type and 
acceptor-type defect). (a) Linear. (b) Semi-log. 
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neutron irradiation, which is consistent with the theoretical 
analysis above. 


IV. NEUTRON IRRADIATION EXPERIMENT 


This paper uses the D-T neutron generator for irradiation, 
which is based on the accelerator neutron source and affiliated 
with the China Academy of Engineering Physics (CAEP), 
Mianyang, China. The experiment is carried out on the K- 
400 accelerator, using accelerated Deuterium (D) beam to ir- 
radiate Tritium (T) target, and D-T reaction occurs, producing 
fast neutrons with energy of about 14.1 MeV. 

To improve the efficiency of neutron irradiation experi- 
ment, a new device layout mode, called stepped layout mode 
is proposed and shown in Fig. 15. Six ABB 5STP 52U5200 
thyristor devices are distributed in three spherical spaces with 
different radii for irradiation, and two thyristors in each radius 


face the neutron source with anode and cathode, respectively. 

According to the neutron calculation report of ITER, the 
neutron flux rate at QPS could be up to 10°n x cm7? x s7! 
with local shielding. In this paper, the average neutron flux 


4s Fig. 12. It can be seen that when the thyristor is subjected 
459 to the same voltage, the forward and reserve leakage current 
40o 1p RM; 1RRm increases under both donor and acceptor defect 
4851 effects. 
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Fig. 15. (Color online) Schematic diagram of neutron irradiation 
experiment. (a) 2D space distribution. (b) 3D space distribution. (c) 
Neutron source. 


rate is taken as 104n x cm7? x s71. Assuming that the duty 


cycle of ITER operating is 0.5, the operating time of a year 
is estimated as 1.58 x 107s, so the neutron flux at QPS could 
reach 1.58 x 10!4n x em~?. In order to obtain more ex- 
treme results, the neutron flux at 20 cm away from the neu- 
tron source is set to 4.3 x 10!!n x em~2, which is 2.7 times 
larger than the actual irradiation environment. Then the out- 
put number of the neutron source can be calculated by the 
following formula. 


NTotal = $ X Amr? = 2.16 x 1015n, (35) 


In the irradiation experiment, the cumulative irradiation 
time was 95373 s, and the final number of neutrons produced 
by the D-T neutron generator is 2.55 x 101°n (uncertainty is 
4%), which meets the experimental requirement. 

In addition, the neutron flux of the devices from near to 
far is 1:1/4:1/9, and the Cross-sectional area of the irradiated 
thyristor is 233 cm?, So the neutron flux and the neutron in- 
jection number of the thyristors at each position could be cal- 
culated. 

As the thyristor stack in QPS is used as the trigger unit of 
commutation circuit, which is shown in Fig. 1, the temper- 
ature of the thyristor is not obvious after flowing through 1 
ms/100 kA pulse current. Hence, the electrical parameters 
that need to be tested when the device temperature is 25 °C 
are taken as the key test data. The tested thyristor number and 
location are shown in Table 1. 


TABLE 1. The thyristor number and location in 14.1 MeV neutron 
irradiation experiment. 
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Fig. 16. (Color online) The change of the thyristor each parame- 
ter before and after neutron irradiation. (a) Change of Iprm. (b) 
Change of Irrm. (c) Change of Urm. (d) Change of Qrr. (e) 
Change of trr. (f) Change ratio. 


s2 0.48 x 10''n/cm? respectively. The forward and reverse 
523 leakage current [pry and IRrm, which represent the block- 
524 ing characteristic of the thyristor. [pray increases to 637 uA, 
55 528uA and 452.54A respectively, and [pray increases to 
se 615.5uA, 458uA and 422.5uA respectively. The on-state 
527 Voltage drop Urm of the irradiated thyristors at different dis- 
ses tances from the neutron source increases from 1.26 V to 4.325 
52 V, 1.84 V and 1.415 V respectively. In addition, the value of 
530 Qrr and t,.,. have a significantly reduce after neutron irradia- 
531 tion shown in Fig. 16d and Fig. 16e, which means the severity 
532 Of the reverse recovery process is reduced and the duration is 


60 on shortened after neutron irradiation. 


Location 20 cm 40 cm 
Neutron Flux |4.3 x 107'n/em?/1.07 x 107'n/em?|0.48 x 107 n/em? 
Number #1 #2 #3 #4 


The change of the thyristor each parameter before and af- 
ter neutron irradiation is shown in Fig. 16. It can be seen 
that the tested electrical parameters of the thyristors have ob- 
vious changes after neutron irradiation. As the thyristors are 
placed 20 cm, 40 cm and 60 cm away from the neutron source 
respectively, the neutron flux thyristors received could be cal- 
culated, which are 4.3 x 101'n/em?, 1.07 x 1011n /cm? and 


#5 #6 s34 IIn terms of impact degree, the blocking characteristic and 


535 reverse recovery characteristic of irradiated thyristor have 
sss been seriously affected, which have hundredfold change 
537 Shown in their key parameters. Meanwhile, the conducting 
sss performance of the most severely irradiated thyristors has also 
53 decreased by about 3/4, which is shown by the rise of the on- 
54 State voltage drop Urm. Moreover, it can be concluded from 
541 Fig. 16 that the higher the irradiation flux received by the 
542 thyristor, the greater impact on the thyristor electrical param- 
543 eters, and the greater the performance damage to the thyristor. 


V. SIMULATION ANALYSIS OF THE INFLUENCE OF 
NEUTRON IRRADIATION ON QPS 


544 
545 


sẹ QPS is a crucial system for not only superconducting mag- 
nets, but also fusion device. It is impractical and risky to 
conduct the irradiation experiments on actual system. With 
the development of simulation technology and PLECS soft- 
ware, the simulation accuracy has been greatly improved, 
which could effectively simulate the real QPS electrical per- 
formance in normal and operating states. Hence, the simula- 
tion method is used to study the impact of neutron irradiation 
on QPS. 

The influence of high-doze neutron irradiation on thyristor 
stack may bring huge risk of breakdown to QPS in fusion de- 
vice. Hence, it is of great significance to analyze the impact 
of neutron irradiation on QPS. By irradiation experiments on 
QPS is risky and unreasonable, so simulation method for pre- 
diction is adopt in this paper. The simplified circuit of QPS in 
s PLECS software and the typical working current waveform 
are shown in Fig. 17 [36]. 
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Fig. 17. (Color online) (a) Simplified simulation circuit. (b) Typical 
current wave of QPS. 


ses According to the results of irradiation experiments, both 
the static parameters IprM, RRM, Urm and dynamic pa- 
rameters Qrr, trr all changed significantly after irradiation 
and may cause the QPS to fail. 

The thyristor stack will be turned off and keep in blocking 
state when QPS in on steady state, in which the capacitor in 
the commutation circuit shall be charged to the preset voltage. 
Therefore, the increase of the [pray may cause the decrease 
of the capacitor voltage. Once the capacitor voltage drops to 
a certain value, the pulse current generated by commutation 
circuit cannot force the current in VCB down to zero, result- 
ing in QPS failure. The larger the Iprm, the short the time 
for the capacitor voltage drops to the threshold. The breaking 
failure waveform is shown in Fig. 18. The capacitor does not 
remain in the reverse blocking state, so the IrRrm has almost 
no obvious impact on QPS. 

The thyristor stack will be turned on and keep on-state con- 
dition only when the commutation circuit generates the re- 
verse pulse current. Therefore, the change of the Urm will 
affect the characteristic of the pulse current, which is indi- 
cated in Fig. 19. It can be seen that the pulse currents before 
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Fig. 18. (Color online) Breaking failure waveform of VCB. 


se and after irradiation are almost the same, which suggested 
sss that the impact of Urm change on QPS is relatively minus- 
ss cule. 
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Fig. 19. (Color online) Comparison diagram of the influence of 


Urm change on the reverse pulse current of QPS. 


Theoretical analysis and experiment verification show that 
sss the reverse recovery characteristics of the thyristor are signif- 
icantly changed. The Q,. and tyr are decreased after irra- 
diation, which means that the reverse recovery speed of the 
irradiated thyristor is accelerated. In QPS, multiple thyris- 
tor devices are connected in series to improve the withstand 
voltage capability. As each thyristor receives different neu- 
594 tron irradiation doze, the change ratio of the reverse recovery 
ss Speed is different. The voltage waveforms of series connected 
sæ thyrsitors are shown Fig. 20. It can be seen that before neu- 
597 tron irradiation, the series voltage sharing effect of the four 
sə thyristors is good. By contrast, after neutron irradiation, the 
sə higher recovery voltage is applied to the thyristor which owns 
eoo the faster reverse recovery speed. This unexpected event will 
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Fig. 20. (Color online) The voltage waveform of the thyristors in 
reverse recovery process before and after neutron irradiation. (a) 
Before neutron irradiation. (b) After neutron radiation. 


cause the device damage due to overvoltage, leading to the 
damage of the entire thyristor stack, ultimately resulting in 
the failure of QPS quench current interruption and transfer 
actions. 


VI. SUMMARY 


This paper studies the influence of neutron irradiation on 
high-power thyristor device in fusion environment. By con- 
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